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Qh' Abstract 

^ ■ In continuum mechanics, the equations of motion for mixtures were derived through 

the use of a variational principle by Bedford and Drumheller in [1978]. Only immis- 
cible mixtures were investigated. We have chosen a different approach. 



c/2 

.y . In this paper, we first write the equations of motion for each constituent of an in- 

£^| viscid miscible mixture of fluids without chemical reactions or diffusion. The theory 

r~| ' is based on Hamilton's extended principle and regards the mixture as a collection 

Q-i< of distinct continua. 

The internal energy is assumed to be a function of densities, entropies and suc- 
cessive spatial gradients of each constituent. Our work leads to the equations of 
motion in an universal thermodynamic form in which interaction terms subject to 
constitutive laws, difficult to interpret physically, do not occur. 
'^ ' For an internal energy function of densities, entropies and spatial gradients, an 

. ■ equation describing the barycentric motion of the constituents is obtained. 

!>■ 

O 

QQ I tained. A form of Clausius-Duhem's inequality which represents the second law of 

^D ■ thermodynamics is deduced. In the particular case of compressible mixtures, the 

equations reproduce the classical results. 

Far from critical conditions, the interfaces between different phases in a mixture 



;_( ' of fluids are layers with strong gradients of density and entropy. The surface tension 



of such interfaces is interpreted. 
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1 Introduction 



The knowledge of equations describing processes occuring in fluid mixtures in 
one or several phases is scientifically and industrially very important ([Barrere 
& Prud'homme, 1973]; [Ishii, 1975]). Many papers concerning theoretical and 
experimental works for mixtures have been produced in continuum mechan- 
ics and physical chemistry. References in [Bedford & Drumheller, 1983] cover 
most of the works in continuum theory. The book by Rowlinson & Swinton 
[1982] gives the statistical methods. 

Truesdell [1957,1965] first derived equations of balance and motion through 
the use of a continuum theory of mixtures. The mixture is then considered 
as a distribution of different continuous media in the same physical space, at 
time t. Each constituent of the mixture is identified with its reference space. 
The thermal and mechanical equations of balance are introduced for each 
constituent. They keep the usual form for a fluid description but they must 
include numerous terms of interaction between the different constituents. The 
study of the average motion, the motion of the barycentre of constituents with 
coefficients equal to their respective densities, are deduced from the sum of 
momenta, energy and entropy of each constituent. The terms including inter- 
action between different constituents of the mixtures are difficult to interpret 
physically. They require constitutive postulates that are difficult to interpret 
experimentally. 

Several assumptions are open to doubt. How must we choose an entropy? 
(Mixture entropy or entropy for each constituent?) What is the form of the 
second law of thermodynamics? Is it symbolized with one or several Clausius- 
Duhem's inequalities? Must we consider a temperature of the mixture or a 
temperature for each constituent? 

MiiUer [1967]; Williams [1973]; Atkin & Craine [1976]; Sampio & Williams 
[1977]; Bowen [1979]; Nunziato & Walsh [1980]; Drew [1983] and many others 
argue about the various points of view. 

In several papers Bedford and Drumheller [1978]; [B & D, 1979]; [B & D, 
1980] presented a variational theory for immiscible mixtures (a continuum 
theory for mixtures whose constituents remain physically separated, such as a 
mixture of immiscible liquids or a fluid containing a distribution of particles, 
droplets, or bubbles). That the constituents remain physically separated has 
several implications for the theory. For instance, the authors use a variational 
principle in which a variation is added to the motion of only one constituent. 
However it is necessary to add to the variations of kinetic and potential en- 
ergies a virtual work due to the forces not represented by a potential. Such 
forces are difficult to take into account in a variational principle deduced from 
the one by Hamilton. It is not possible to use such a method in a more general 
model. 

In this paper a new mathematical systematic method that leads to the equa- 
tions of motion and energy for miscible fluid mixtures without chemical reac- 



tion and diffusion is propounded. The presented work is essentially different 
from the one by Bedford and Drumheller. 

The mixture is represented by several distinct continuous media that occupy 
the same physical space at time t. The novelty consists in the fact that the 
used variational principle is applied to a lagrangian representation associated 
with a reference space for each component. The proposed method uses the 
knowledge of the internal energy of the mixture at every point of the physical 
space. This is not the case in [B & D]. 

For classical fluids, the method corresponds to Hamilton's principle in which 
one makes a variation of the reference position of the particle. It was proposed 
by Gouin [1987]. 

The inferred equations are written in a different form that facilitates the study 
of first integrals for conservative motions ( [Casal, 1966]; [Gouin, 1981]; [Casal 
& Gouin, 19856, 1988a, 1989]). 

For mixtures, the principle appears far more convenient than the method of 
only taking the variation of the average motion into consideration: it separately 
tests each component and leads to a thermodynamic form of the equation of 
motion for each of them. 

The internal energy is assumed to be a function of different densities or con- 
centrations of the mixture. Because of the form of the principle we use, one 
must consider an entropy for each component in the mixture. The total en- 
tropy of the system is the sum of the partial entropies. 

To consider areas where strong gradients of density occur - for example shocks 
or capillary layers - the internal energy is chosen as a function of the successive 
derivatives of densities and entropies. 

Two cases of conservative motions for mixtures are considered: isentropic mo- 
tion and isothermal motion. These are only limiting mathematical cases for 
testing the method. The physical motion is intermediate, between these ex- 
tremes. 

An important medium is the one of thermocapillary mixtures. It corresponds to 
an energy taking into account the gradients of densities and entropies. Then, 
one obtains an equation of the barycentric motion of constituents. For a com- 
pressible mixture, the equations of motion involve a hydrostatic stress. 
In the same way as for thermocapillary fluids, an additional term that has 
the physical dimension of a heat flux may be added to the energy equation 
([Eglit, 1965]; [Berdichevskii, 1966]; [Casal & Gouin, 1985a]; [Dunn & Serrin, 
1985]; [Gatignol & Seppecher, 1986]; [Seppecher, 1987]). A generalization of 
the Clausius-Duhem inequality is obtained for non-conservative mixtures. The 
inequality is deduced from the heat conduction inequality by introducing an 
irreversible stress tensor associated with a general dissipation function. The 
theory flts with the second law of thermodynamics. 

From these results one generalizes equilibrium conditions obtained through 
interfaces in compressible fluid mixtures [Rocard, 1967]. An interpretation of 
capillarity stresses is deduced (rule of Antonov) [Emschwiller, 1964], and gen- 
eralization of results by Vignes-Adler & Brenner [1985] may be investigated. 



2 The motion of a continuous medium 

2.1 Motion of a fluid with only one constituent 

Recall that the motion of a continuous medium can be represented by a sur- 
jective differentiable mapping: 

z^ x=m{z), 

where z = {t, x) belongs to W , an open set in the time-space occupied by 
the fluid between time ti and time ^2- The position of a particle in a reference 
space Vo is denoted by X; its position at time t in Vt is denoted by x [Serrin, 
1959]. Hamilton's principle - variational form of the principle of virtual powers 
- allows us to investigate the equation of motion. 

The variations of motion of particles are deduced from: 

X=*(cc,t;/3) 

for which /? is a parameter defined in a neighborhood of zero; it is associated 
with a family of virtual motions of the fiuid. The real motion corresponds to 
/3 = 0. 

Virtual material displacements associated with any variation of real motion 
can be written [Gouin, 1987]: 



«^f 



/3=0 



Such a variation is dual with Serrin's [S, 1959]. This has been studied in the 
case of compressible perfect fiuids [Gouin, 1978] and corresponds to the natural 
variation of the motion in a lagrangian representation. 

The lagrangian of the fiuid is: 

where V denotes the velocity vector of the particles , Vt the extraneous force 
potential defined inVt , a the specific internal energy, p the density and * the 
transposition in Vf. 

Between time ti and time t2, the hamiltonian action is written: 

/ Ldvdt. (2) 

ii JVt 



The use of virtual displacements (1) provides the equations of motion of com- 
pressible perfect fluids in a thermodynamic form. For fluids of grade n, the 
calculus is swift. It generalizes the equation (29,8) of fluids in [Serrin, 1959] 
that may be written in a form that is independent of their complexity [Gouin, 
1987]: 

r = ^ grad s — grad(/i + Q) , 

where F denotes the acceleration, s the entropy, 9 a temperature and h a 
specific enthalpy. 



2.2 Motion of a fluid mixture 

For the sake of clarity, we study a mixture of two fluids. The method can be 
immediately extended to any number of constituents. No assumption has to 
be made about their composition or their miscibility. 

The motion of a two- fluid continuum can be represented using two surjective 
differentiable mappings: 

z^Xi = Mi{z) and z^ X2 = Mal^), (3) 

(Subscripts 1 and 2 are associated with each constituent). 




Fig. 1. 

Xi and X2 denote the positions of each constituent in reference spaces V^i 
and Vo2- The variations of motion of particles are deduced from: 

Xi = *i(a;,t;/5i) and X^ = ^2ix,t; [3^). (4) 

This is a generalization of the one in paragraph 2.1 where /3^ and (^2 ^^^ defined 
in a neighborhood of zero; they are associated with a two-parameter family 



of virtual motions of the mixture. The real motion corresponds to /^^ = and 
j32 = 0; the associated virtual displacements can be written: 






dp. 



and 82X2 



f3,=0 -9/3, 



(5) 



They generalize relation (1) given for a fluid with one constituent. The prin- 
ciple tests separately each constituent of the mixture and so is very efficient. 

Remark. - It follows from the classical representation of a two-constituent 
mixture, that formula (4) can be written 

x=ifi-,{Xi,t;P,) and a; = (^2(^2, t; /^s)- (6) 

A variation of the motion like Serrin's (p. 145, [S, 1959]) cannot be extended 
to any two-parameter family of virtual displacements: as a matter of fact, 
any field of virtual displacements associated with (f, must be such that 5x = 

' Xi, t; 0) and there is no reason for being identifiable with a vector field 



dp 



dip2 



written in the form 5x. = — -— (X2, t; 0). 
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In the case of compressible fiuids, writing (1) for the virtual displacement 
associated with a one-parameter family of virtual motion of the fiuid is deduced 
from the variation defined by Serrin by the way of vector spaces isomorphism 
[Gouin, 1976]. For mixtures, by using (5), a virtual displacement of miscible 
fluid mixtures can be considered. 

The lagrangian of the mixture is: 

L = -^Pi'^l^i + 2P2V2V2 -e - pi^i - p2Vt2. 

Here, Vi et V2 denote the velocity vectors of each constituent, p, and p2 are 
the densities, Vti and VL2 are the extraneous force potentials depending only 
on z= (t, x) and e is the internal energy per unit volume. 
The expansion of the lagrangian is general. In fact dissipative phenomena 
imply that Vi is almost equal to V2. Because of the interaction of the con- 
stituents, e does not always divide into energies related to each constituent of 
the mixture, as it does in the case of simple mixtures of fluids [Miiller, 1968]. 
More generally, the lagrangian may be written 

L = ^PiVtVi + ^P2V;V2 - 11 

where r] denotes the volumetric potential energy of the mixture. In our par- 
ticular case, T] is the sum of internal energy and extraneous force potentials. 
In fact, the following calculus would be similar (but even simpler). We use the 



first form of the lagrangian in order to obtain classical expressions used for 
fluids. 

Hamilton's action between time ti and time t2 can be written as in (2). The 
mixture is assumed not to be chemically reacting. The conservation of masses 
requires that: 

p,detF, = p„{Xi) (7I) 

for the densities of each constituent, with subscript i belonging to {1,2}. At 
fixed t, the jacobian mapping associated with each Mj is denoted by Ff, p^^ is 
the reference specific mass in T?oi. 
In different iable cases, equation (7^) is equivalent to: 

^ + divp,V, = 0. (7^) 

The volumetric internal energy e is given by the thermodynamical behavior 
of the mixture and will appear in the equations of motion. Each constituent 
has a specific mass; similarly, two specific entropies Si and S2 are supposed to 
be associated with constituents 1 and 2. 

For a specific internal energy depending on derivatives of gradients of densities 
and entropies up to the n-th order: 

e = e(si, grad si, ..., (grad)''"^si; S2, grad S2, ..., (grad)""^S2; 

Pi,gradpi, ..., (grad)""Vi;P2>gradp2, •••, (grad)""V2)) (8) 

(we say, the mixture is of grade n), the calculus is similar to than in ([Gouin, 
1987]; [Casal & Gouin, 1988a]). The orders of gradients of pi, p2, si, S2 can 
be chosen differently, but the exposition would be less clear. Let us note that 
the volumetric potential energy can be written: 

?7 = r7(si, grad Si,...,(grad)""^si;s2, grad S2,---,(grad)""^S2; 
Pi, grad pi, ..., (grad)"' Vi; P2> grad P2, •••, (grad)"~V2) ^) 

with z = {t, x). 

2. 3 Examples of motions of a fluid m,ixture 



Let us consider two examples of Hamilton's principle. Obviously, they are two 
limit cases deduced from a variational principle and incompatible with any 
dissipative process. 

First case: Isentropic motions. 



We assume that the specific entropy of each particle is constant: 

^ = or — + —V- = 
dt ' dt dx ^ 

(Lie's derivative with respect to the velocity field of constituent i is equal 
to zero. The average motion of the mixture does not have Lie's derivative, so 
henceforth, we can expect each constituent of the mixture to possess entropy). 
The equation 



^oi \ -^i 



(9) 



defines an isentropic motion. 

Motion equations are deduced from Hamilton's principle. The two variations 

of hamiltonian action are: 



Sitt = a'{f3j) 



/3,=0 



For i G {1, 2}, we deduce: 



Sia= I I I (^'^*i^i - £,p, - ^ij ^iPi + AV*(5iVj - e^s.ksi 



5iSi^ + ... + £., 



S7 t7 



5jSj 

>,71>...>7„-1 ' ',71,...,7n-l 



^Pt,^ ^r^l,-l '^»'71,...,7„-l '^*.7l,...,7„-l 



dvdt. 



Subscript 7 corresponds to spatial derivatives associated with gradient terms; 
as usually, summation is made on repeated subscripts 7. 



The quantity 



p=i 



defines the specific enthalpy of the constituent i of the mixture, and 



(IQi) 



1 



Pi 



e 



E(-ir(^,«. 



p=i 



(10^ 



defines the temperature of the constituent i of the mixture. 

The motion of constituent i of the mixture is expressed by the equation 

Tj = 6'igradsj -grad(/ij + fij). (11) 

Appendix 1 gives details of the calculation. 

The fact that an entropy is introduced for each component implies the exis- 
tence of a temperature for each fluid constituent. This is in accordance with 



[Green and Naghdi, 1965]. In this first limit case, the variational principle 
does not lead to the equality of temperatures. The behavior of the two com- 
ponents corresponds to their thermodynamical independence (no exchange of 
heat between particles). The two components are only bounded by effects of 
compressibility. Such a model can be used in practice only if one includes 
the fact that some heat conduction exists and adjusts the values of the two 
temperatures. So, one verifies that equations obtained with the variational 
principle must be completed. 

Second case: Only the whole entropy of the mixture is conserved. 

We have: 

/ iPiSi + P2S2)dv = K, (121) 

JVt 

where K is constant (independent of the time). 

Let us use the method applied in [Casal & Gouin, 1988 &] for compressible 
fluids. The variation of the entropy of each component is assumed to be the 
sum of a variation associated with the virtual motion and another one, 6oiSi 
related to the particle: 

dXi 

where SoiSi and 5iXi are independent and 5oi must verify, with respect to 
component i of the mixture: 

rt-z 



^oi / / (PiSi + p2S2)dv = 0. (12^ 

Jti JVt 



Glassical methods of variational calculus provide the variation of the hamilto- 
nian action: the variations 6iXi yield a motion equation of the form of (11). 
Relation (12^) implies: there exists a constant Lagrange multiplier 9o, constant 
in space and in time, such that: 



rt2 



ti JVt ^ n " 



<^ojSi,^j,...,^^_^ + PiOo^oiSi) dvdt. 



'*i,71,...,7„-l -"'■'',71,.. .,7. 

Let us point out that adiabatic condition (12^) implies that even if the mix- 
ture is moving, 9o remains constant in space and time since 6o is a Lagrange 
multiplier associated with an integral constraint. 

Considering variations 5oiSi vanishing on the boundary of Vt ; an integration 
by parts yields: 

Soia= / pi{6o-0i)SoiSidvdt. 



ti JVt 



The principle: "for any 5oiSi vanishing on the boundary of P^ , ^oja = " leads 
to the equation for the temperature: 

e^ = e,, (13) 

i.e., all the constituents have the same temperature 60, constant in all the 
flow: the flow is isothermal. 

This variational principle can be compared with the static one expressing 
that the equilibrium of a compressible fluid kept in a fixed adiabatic reservoir 
without exchange of mass and entropy with outside is obtained by writing the 
uniformity of temperature and pressure. It can be interpreted as a limit case 
of a mixture with an infinite heat conductivity. 
Equation (11) is replaced by 

r, + grad((7, + Q,) = 0, (14) 

where gi = hi — OiSi is the free specific enthalpy (or chemical potential) of 
constituent i. 

When the fluid is compressible, Eqs (11) or (14) yield first integrals and theo- 
rems of circulation [Serrin, 1959]; [Casal, 1966]; [Gouin, 1981]; [Casal & Gouin, 
1985 &, 1989]. These properties can easily be extended to each constituent of 
an inviscid mixture. 
Each constituent of the mixture in an isothermal equilibrium state satisfies 

G, = Q (15) 

where Gi = Qi + Qi and Cj is a constant. 

An interface in a two-phase mixture is generally modelled by a surface without 
thickness. Far from critical conditions, this layer is of molecular size; density 
and entropy gradients are very large. A continuous model deals with such 
areas by using energy in the form (8) extending forms given in ([Rocard, 
1967]; [Cahn & Hilliard, 1959]) for compressible fluids. 

Equation (15) is satisfied for each constituent of a mixture in an isothermal 
equilibrium state. The fact that Cj values are the same in each bulk is also a 
consequence of Eq.(15). It follows that with extraneous forces being neglected 
in a two-phase mixture, if the two different phases are in isothermal equilib- 
rium, then specific free enthalpies (or chemical potentials) of each bulk are 
equal. This result generalizes the classical law to mixtures of grade n: the 
different constituents keep their own chemical potential through a change of 
phases [Rocard, 1967]. 
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3 Thermocapillary mixtures 



The previous equations yield the equations of motion for each of the two fluids 
constituting the mixture. Usually, the chemical potentials hi and Qi are rarely 
considered except for equilibrium; Eqs. (11) and (14) make use of them in the 
study of the motion. 

In order to compare results given by the variational technics and the usual 
approach of mixture motions, we can try to study the mean motion and diffu- 
sive motions for two components. But how do we define these motions? In this 
paper, we consider the mean motion using an approach which is independent 
of the number of constituents. 

For the mean motion, p = J2iPi and p T = J2i Pi^i (using T = Y^i'^i has no 
physical meaning). 

On the contrary, p T has the physical dimension of a volumetric force and 
is an additive quantity in the same way as body forces or the stress tensor 
divergences. So, in this first article, we have not attempted to investigate the 
difficult question of the relative motion associated with " Fi — r2 " and hence 
to see how diffusion behaves in miscible mixtures. 



3.1 Conservative motions for mixtures of thermocapillary fluids 

In the case n = 2, we write the equations of motion by using the stress tensor. 
The mixture is said to be a thermocapillary mixture if the internal energy e is 
written as: 

e = £(si,gradsi,S2,grads2,Pi,gradpi,p2,gradp2)- (16) 

The expressions 

2 2 

i=l i=l 

introduce the total density and acceleration of the mixture. With denoting by 
a the specific energy of the mixture, we have 



and, 



with 



pa = e, 



2 'T^ 1 1 

da = y^^—^dp.: + —Qi dsi + -^i^dpi ^ + -"^i^dsi^^, 
f^iPPi P P P 

Pi 



Four new vectors $j and ^j {i G {1, 2}) are introduced by the theory; it is an 
extension of the calculus written in [Casal & Gouin, 1988a]. 
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For a compressible fluid (with only one constituent), V and O are respectively 
the pressure and the Kelvin temperature. The fluid is isotropic, so the internal 
energy is of the form 

^ ^[Pii "^J? Pij^ Xiji lij) 

with (3ij = giadpi.giadpj, Xij = gradp^.grads^, 7^^- = grads^.grads^ (i e 

{1,2}). 

It follows that: 



$ j = ^ Cij grad pj + Eij grad Sj 

2 
*j = XI ^ij Srad Sj + Eij grad Pj, 

with 

Cij = (1 + Sij)e^i3^., Dij = (1 + 5ij)£,7,^., £^ij = £,Xy, 

where (5^^ is the Kronecker symbol. The simplest model is for Cij, Dij and Eij 
constant. 

Assuming that f^i = ^2 = f2 (this is the case for gravity forces), we obtain a 
general formulation of the equation of motion for thermocapillary mixtures: 

pF = diva — pgradn (17) 

or 

pT^ = ay.,^^ - pQ^^. 

The sum of ai and o"2 obtained by successive integrations by parts is conven- 
tionally called the generalization of the stress tensor. This is simply for the 
convenience of writing (in fact, the sum of stresses associated with different 
media is a non-sense). 

Let 

2 

i=l 

where 

CTiu^ = -piSuy - ^iuPi.^ - "^iuSi^^, (18) 

with 

Pi = Vi - pidiv^i. 

The tensorial natures of entropy and density are different; that is the reason 
why only div # appears in the pressure term. The calculus for this case appears 
in appendix 2. In appendix 4, we deduce Antonov's rule related to the surface 
tensions of an interface between two immiscible fluids from the expression of 
the stress tensor. 

For mixtures with internal capillarity ([Gatignol & Seppecher, 1986]; [Seppecher, 
1987]), only $j, i G {1, 2}, are not equal to zero. In the case of classical com- 
pressible mixtures, terms of gradients are null and $i and ^j are zero; a^^ is 
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written: 



-V5,^ 



with V = J2'i=i ^i defining the pressure of the mixture. 

Let us remark that Vi with i G {1,2} does not correspond to the classic notion 
of partial pressures in ideal mixtures. This result fits with I. Miiller's paper 
([1968], p. 37): according to Miiller terminology, ii s = pa = piai + P2<^2 with 
ai = ai{pi,Si), such an internal energy corresponds to a simple mixture of 
fluids. If IIj = pfaip, (pj. Si), then 

Pi = ni + ^i^(«i-a2) 



and 



P2 = n2 + ^^(a2-ai) 
P 



with assuming that for every point in the mixture all the constituents have the 
same temperature 6i = 9q, then, P = Hi + 112 and IIj represents the partial 
pressure of constituent i [Bruhat 1968]. 

Example of binary system with complete miscibility: liquid-vapor equilibrium 

At a given temperature, the numerous constitutive equations represent V as 
a function of densities of the constituents ([Soave, 1972]; [Simonet & Behar, 
1976]; [Peng & Robinson, 1976]). These equations, very useful in petroleum 
industry, derive from subtle calculus on mixtures of Van der Waals fluids. 
Taking into account Eq. (15), they allow the study of a liquid bulk and a 
vapour bulk for two fluids that can be completely mixed. 
A plane interface lies between two phases of a two-fluid mixture in isothermal 
equilibrium. This is an unidimensional problem with reference to the direction 
perpendicular to the interface. Body forces being neglected, Eq. (17) can be 
written: 

dz 
Vectors $j and ^j are zero in the bulks in which Vi = Vy = Vq (/ and v refer 
to the bulks with the same pressure Pq). 

For equilibrium, eight variables p^, p^^, su, Siv with i G {1,2} determine the 
properties of the two bulks. They are subject to six relations: 

9ii = giv, 
921 = g2v, 

dll = dlv = 021 = d2v, 

and 
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This fits with Gibbs law of phases that gives the system as " divariant" [Rocard, 
1967]. Two additional conditions, On = 9iy = Oq {i E {1, 2}) and Vi = Vy = Vq, 
set an equilibrium of the mixture. 



3.2 Equation of motion and equation of energy for a thermocapillary mixture 



Conservative motions 

Some properties of mixtures are directly deduced from the equation (17) of 
motion. 



Let us define 



Mj = Ti - 9i grad Si + grad(/ii + Qi), 

5. = ^ + AdivV„ 
at 

Q - 0^^' 

Ei = —^ + div [(ci - ai)Yi] - div Uj - p^ 



dt ' ""■ ''^'' "'^ "'' "'■ "' " dt ' 

with 

e^ = ^ aV2 + e, + p^(]„ £, = ^ 
z p 



and 



dPi ^ dsi ^ 
Ui = -p $i + — ^ *i 
dt dt 

MPi ^ ^ ^ ^y. g^^^ ds^^ds^^ ^V- I 

Y dt dt dx * dt dt dx y 

Theorem 1 For a thermocapillary internal energy and for any motion of the 
mixture, the relation 

Y.E,- pJs/l*Yi - (\yI + K + n^ B,-S, = (19) 

is identically satisfied. 

(Proof in appendix 3). Limit cases (the isentropic and isothermal motion of 
Sec. 2) yield the two corollaries: 

Corollary 2 - Any isentropic motion of an thermocapillary mixture satisfies 
the equation of balance of energy: 

j:^ + div[(e, - a,) V,] - div U, - p,^ = 0. (20) 
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with fi = - Pi^l + V'j + Pi^i and (p^ = e^ — PiOiSi corresponding to the specific 



This results from the simultaneity of equations Si = 0, Bi = and Mj = 0. 

Corollary 3 - For any conservative motion of a thermocapillary perfectly heat 
conducting mixture, Eq. (13) is equivalent to: 

E ^ + div[(/. - a.)V,] - divU. - P.^ = (21) 

1 

2 

free energy of constituent i at temperature 6q. 

Equation (21) is called the equation of balance of free energy. 

Dissipative motions 

If we try to introduce dissipative phenomena, while keeping the form of the 
equations, we may only add an irreversible stress tensor an for each con- 
stituent. So, for each fluid, {i G {1,2}), supplementary energy flux vectors 
div(o'/jVj) and dissipative forces tr{ajiA.i) have additive properties (A, is 
the deformation tensor of constituent i) and are introduced and summed. 
The only assumption for functions of dissipation, I]i=i tr^an Aj) ^ 0, leads 
to a Clausius-Duhem inequality. So, the model fits with the second law of 
thermodynamics independently of any constitutive behavior for an and <t/2 
(Chapman and Cowling [1970]). 

Equation of motion (17) can be written: 

pT = div(o" + o"/) — p gradf2 

(as in Eq.(17), we are led to formally writing ai = an + cr/2). 

Afterwards, the heat flux vector will be denoted by q and heat supply by r. It 
follows. 

Corollary 4 - For any motion of a thermocapillary mixture, the equation of 
energy 

^ de- dil 

E ^ + div[(e, -a,- a,,) V,] - div U, - p,-^ + div g - r = (22) 

is equivalent to the equation of entropy.- 

2 dSi 

E Pi^i~J7' ~ tr{aiiAi) + div g - r = 0. 
Assuming that the work of dissipative forces satisfies Z^Li i^^(c"/j^i) = 0, we 
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obtain a generalized form of Planck's inequality: 

2 dSi 

At any point of the fluid, each constituent is supposed to have the same tem- 
perature (not necessarily uniform). This requires that the time for relaxation 
of molecules is short in comparison with the time characterizing the flow. 
Then, 

By use of Fourier's law in the form: 

g.grad^ ^ 0, 
we deduce a generalization of Clausius-Duhem's inequality: 

This Clausius-Duhem's inequality is a form of the second law of thermody- 
namics. Let us note that some authors separate q and r in qi and rf, this does 
not seem necessary. 



Appendix 1 

Equations of motion for mixtures of perfect fluids 

Isentropic motion 

Let us denote by ^j the variation diXi, we deduce from the two relations in 
lagrangian variables [Gouin, 1987]: 

<5.V, = -F.^, S,p, = p, div,, ^, + A|^ i^ (23) 

and 

Os ■ 

dXi 

d 
where — — is the linear form associated with the gradient and divoj is the 

divergence operator on Voi. 

Variations of entropies satisfy relation (9). With assuming that terms on the 

edge of W are zero, we obtain: 
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Sia= I I I [2^i^i - ^,ft - ^i) ^iPi + Pi^i^i^i - ^,s,SiSi 







"n-l 








S(-i)''(^'..„, ,).. 

"n-l 






— 


_p=l 


(5iSj > dvdt, or 


Sitt 


jti . 


/ Poi 




^' ' ' dt ^dxr'\ 



dvnidt 



with Rj denotes -V*Vj — hi — Vti. 



But, 



So, 



n „ art 

diVoi(PoiRi^i) = PoiRidiv^i |. + Ri^ li + Poi^ ^i- 



Sitt 



t2 



ti "'X', 



Po 



s^' . « (v-F.) - e, ^'' 



dXi dt 



dX,, 



$,i dvoidt, 



(terms given by integration on the edge of Voi are zero). 
Finally, for each constituent. 



dt^ '^^~ dXi 'dX,' 



and by using of the relation 

d 



,aV,: 



('V*F-') = r*F- + V* 



we obtain 



d 1 (9R- ds- 



or 



' ' ' dx.,'2'' "'' dXi ' "'dXi 

This leads to the vectorial form: 

Fj = 6i grad Sj — grad(/ii + Vti 



dsi d 
'' dx dx 



ihi + rii). 



:iii 



Appendix 2 
Equation of motion of thermocapillary mixtures 



Equation (11) implies 



J2Pi^iy + Pi^in = E^A^M - (^A,J,.Si,^ -Pi(e,ft),^ +Pi[(£,ft,J,J„ 



j=l 



i=l 
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By noting that 



^,7 = H^A^M + (^A.J,.Si,.7 +^,P^Pi,^+^.P^,.Pi,.^^ 



i=l 



we obtain 



or 



i=l 

I ^ 

U=i 

2 2 

i=l 1=1 '"' 

In the case VLi = 0.2 = ^, this agrees with Eq. (17) when pa = e. 

Appendix 3 



Proof of relation (19) 

Relation (10^) yields 

Ei = -^ + div 
ot 



P. ( 2 V' + /i^ + fi. ) V. 



— div 



dpi ds 



dt 



*-+sr* 



Pi 



dVLi 
'dt 



On the left hand side of Eq.(19), the extraneous force potential and inertia 
terms cancel out. 

With ei,hi,^i and ^j having been replaced by their respective values, it 
remains to be proved that the terms from the internal energy also cancel out. 
By using the convention that terms subscripted by i are summed on {1,2}, 
the internal energy involves: 

a) in Ef -^ + dw{p,hiVi) - div (^^i + ^*i j , 

b) in piM*Vi : -p^Oi grad* Sj Vj + p^ grad* hi Vj , 

c) in Q V^ + h, + n,>jB,:hJ^ + p, div V, j , 



d) in 5 : pfi. (^ 



This leaves the following in the first term of (19): 



I + div(,...V,)-div(f *, + §*, 



-Pi^i [-g^j - Pi grad*/ii ^i~^i[~^ + Pi ^i^ ^i ) ' 



or: 



9^ (u , A- if.\^Pi ( Q ^ A- ,T, ^^^^ <9gradp^ ggrads^ 

--(/.. + div$.)— -(pA + div *.)--$. ^^ *^^^- 

Definitions (16) of e and (10) of hi and 9i show immediately that this expres- 
sion is identically zero. 



Appendix 4 



Interpretation of capillary tensions associated with two immiscible 
fluids 



We consider an interface between two immiscible fluids. The assumption that 
two fluids remain absolutely separated is not reflected in physical experiments. 
The two fluids interpenetrate each other in a thin layer (thickness about few 
molecular diameters). Energetic reasons justify the assumption [Brin, 1956]. 
The density of each constituent regularly changes between the value in bulk 
1 and the one in bulk 2, with a profile in accordance with the one of a liquid- 
vapour interface; each density vanishes asymptotically in the complementary 
phase (see Fig. 2). 

In the interfacial layer, density gradients are important. In the isothermal case, 
the zone of mixture of two isotropic fluids can be represented with an internal 
energy of the form: 

e = pa{p^,p2,f3^,f32,-f) 
with 

/3j = (gradpj)^ and 7 = grad p^ grad p2 • 

The stress tensor reduced from Eq. (18) is: 



a = ^ [-pPi«,p^ + Pi div(Cigradpi + Dgradp^) 

i¥=j 

—Ci grad Pi gr ad* Pi — D grad p, grad*Pj , 



Id 
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with Ci = '2pa^i^, and D = pa^^. Let us denote 

n = XI PPi'^,P, - Pi div(Ci grad p^ + D grad pj + d (grad pj^ + DgradV^ gradp^- . 

In each bulk, gradient terms are zero and 11 takes the values: 

n, = pja,Jp„0,0,0,0) and U^ = p^a^jO, p^,0,0,0) 
that correspond to the pressure in each phase of fluid. 



Bulkl 




Thickness of the interface 



Fig. 2. Profiles of densities through the interface 

When the extraneous force potential are neglected, 

diver = 0, 

which is the equation of the equilibrium state of the mixture. 

For a flat interface normal to gradp^ and gradp2, the eigenvalues of stress 
tensor a are 

Ai = -11 + Ci(gradpJ^ + C2(gradp2)^ + 2L'grad*pigradp2, 

(associated with the plane of interface), and 

(associated with the direction normal to the plane of interface). 

In a system of coordinates suitable for the interface {i.e. an orthonormal sys- 
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tern whose third direction is z, the stress tensor is written: 



a 



Ai 











Ai 











A2 



The equation of balance momentum in the planar interface implies: 

A2 = -Ho, 
where Ho denotes the pressure that must be common to the two bulks. 
The line force per unit of length on the edge of the interface is (see Fig. 3): 

J-Z2 

JF = / Xidz = 

J Z\ 

-Iio{z2 - zi)+ / {Ci(gradpi)^ + C2(gradp2)^ + 2Dgrad*Pigradp2}c^^, 

where Z2 — zi corresponds to the physical size of the interface. 
Let us note that, due to the thickness of the interface, — no(-22— -Zi) is negligible. 
Conditions on the edge are not the ones of first gradient fluid. Nevertheless, 
the problem concerns a plane interface and the mean radius of curvature Rm 
is such that Rm = 00. Therefore, the constraint condition on the edge is in 
the form n*a where n* is normal to the edge of the interface [Casal, 1972]. Let 

m 




Fig. 3. Interpretation of the mixture surface tension 
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us note: 

Hi = {Cj(gradpj) } dz et i^i.2 = / 2Dgia,d*p-^^gra,dp2dz. 

J Z\ J Zl 

(The integrals converge in the two bulks). The line force per unit of length on 
the edge is: 

H = Hi + H2 + Hi2- 

Here, H represents the surface tension of the flat interface in its equilibrium 
state expressed as a function of the surface tensions Hi and if 2 of each bulk 
and a supplementary term ifi.2 associated with the effect of the mixture. This 
is the expression of the surface tension of an interface between two immiscible 
fluids known as Antonov's rule [Emschwiller, 1964]. 
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